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We calculate the critical temperature of the Ising model on a set of graphs representing a concate- 
nated three-bit error-correction code. The graphs are derived from the stabilizer formalism used in 
quantum error correction. The stabilizer for a subspace is defined as the group of Pauli operators 
whose eigenvalues are +1 on the subspace. The group can be generated by a subset of operators 
in the stabilizer, and the choice of generators determines the structure of the graph. The Wolff 
algorithm, together with the histogram method and finite-size scaling, is used to calculate both the 
critical temperature and the critical exponents of each structure. The simulations show that the 
choice of stabilizer generators, both the number and the geometry, has a large effect on the critical 
temperature. 
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I. INTRODUCTION 

A bit of information can be stored in any physical sys- 
tem with two distinct states. For the physical system to 
be a reliable memory, the states of the system must be 
robust against external fluctuations. The classic example 
is a ferromagnet. Below a critical temperature, the size 
of the average magnetization is robust against changes 
in external magnetic field and temperature-driven spin 
fluctuations. 

A quantum bit (qubit) of information can be stored 
in any physical system with two orthogonal quantum 
states. A goal of quantum information is to engineer 
a system that can reliably store the state of the qubit 
in the presence of environment-induced fluctuations. A 
number of approaches have been proposed from quantum 
error-correction [l| to passive protection of the informa- 
tion through symmetries Q or energetics [1, d, H, H, 0] • 

The possibility of engineering the quantum equiva- 
lent of the magnetic hard drive is quite attractive. The 
premise is that a macroscopic number of qubits with 
multi-qubit interactions could create a single stable qubit 
memory. It is widely suspected that Kitaev's toric code 
on a four dimensional lattice would achieve this task 
(HQ. In lower-dimensions, the answer is unclear. Bravyi 
and Terhal have recently shown that for interactions 
based on stabilizer codes, there is no two-dimensional 
self-correcting quantum memory Q • They make the rea- 
sonable physical assumption that the number of qubits 
involved in the interactions does not grow with the size of 
the lattice. In the case of self-correcting memories based 
on concatenated codes the number of qubits involved in 
each interaction does grow with the lattice size pj ■ 

Here we consider the classical concatenated triple- 



modular redundancy code in the formalism of quantum 
stabilizers. The standard choice of generators for this 
code leads to interactions that grow with the system 
size. Choosing a different set of generators yields in- 
teractions that are only between two-bits and equiva- 
lent to an Ising model. In this paper, Monte Carlo 
simulations are used to study the critical behavior of 
the error-correction inspired structures shown in Fig. [T] 
within the framework of a ferromagnetic Ising model. 
The calculations use the Wolff algorithm [lCj together 
with the histogram method [ll], [l2| and finite-size scal- 
ing [l3j]. These high resolution Monte Carlo techniques 
have been utilized successfully to study the critical phe- 
nomena of many different model Hamiltonian systems 
such as the 3D Ising ferromagnet Heisenberg lat- 
tice [Hi, XY models QJ, [T3|, dilute Ising magnet [II], 
Potts models [l9j], and Sier pins ki fractals of dimensions, 
d, between one and two [20l.l2ll. |22| and between two and 
three [H[ . We study how the choice of generators changes 
the Hamiltonian and affects the magnitude of the critical 
temperature. Critical behavior is characterized by the 
set of critical exponents (a, [3, 7) and Wolff dynamical 
critical exponents are calculated for each structure. 

The results presented here show that structures with 
low dimensionality and two-body interactions preserve 
one bit of information. It suggests two new directions for 
examining self-correcting quantum memories: 1) choos- 
ing non-standard stabilizer generators to minimize the 
multi-qubit interactions and 2) to examine stabilizer 
codes on fractional dimensional geometries. 

II. MODEL AND METHODS 
A. Stabilizer to Structure 
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We consider the familiar classical code of triple modu- 
lar redundancy. Each bit, x, is encoded into a logical bit, 
x L , consisting of three bits of equal value, i.e. L = 000 
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and 1 L = 111. If an error occurs on a single bit, majority 
vote can be used to determine the value of the logical bit; 
two errors would not be corrected. One way to protect 
against higher errors is to concatenate the code recur- 
sively. At each level of concatenation k, the logical bit 
consists of three bits of level k-1, e.g. a at level k = 2 is 
defined as fc = 2 = o fe=1 O fe=1 O fe=1 = 000000000. Correc- 
tion works by majority vote at the lowest level first and 
then working up. Each level of concatenation k can al- 
ways correct a maximum of 2 k — 1 errors on the physical 
bits. 

The basic idea of the stabilizer formalism is that a 
quantum state or subspace can be described b y th e op- 



erators that have +1 eigenvalue on that space [24j. The 
stabilizer formalism is particularly useful for describing 
quantum error correcting codes. Classical error correct- 
ing codes represent a subset of quantum error correcting 
codes that only protect against classical bit-flip errors 
and not the phase errors. The triple modular redun- 
dancy code is a textbook example for introducing the 
idea of stabilizer error correcting codes [23 |. Below we 
review the case of this concatenated code and show how 
it translates directly onto the Ising model. 

Following standard quantum computation notation, 
the states of the j th spin are represented as |0) ■ and |1) ■, 
Zj is the Pauli-z operator on the j th spin (Zj \0} j = \0) i 
and Zj = -|1} 7 ), and X, 
tor on the j th 



J 3 <~>1 '"'J 

, .ij is the Pauli-x opera- 

spin (flips the j th bit). For level- 1 en- 
coded bits |0} fc=1 = \0) 1 |0) 2 |0) 3 = |000) and |l) fc=1 = 
|1) 2 |1) 3 = | 111) . The encoded Z operator is de- 
fined as Z k = Z\~ x Z\~ x Zl~ x and for level 1 is Z k=1 = 
Z\Z 2 Z%. 

For stabilizer codes, the stabilizer is defined as all the 
products of Pauli operators that act trivially on the code 
space. For the three-bit code, the stabilizer consists of 
four operators S k=1 — {I, Z\Z 2 , Z1Z3, Z2Z3} where I is 
the identity. The stabilizer can also be defined by the 
generators of the group, S k=1 =< Z\Z 2 , Z 2 Z 3 >. 

In the case of level- 1 triple modular redundancy, the 
freedom in the minimal generators is trivial S k=1 =< 
Z\Z 2 , Z2Z3 >=< Z1Z2, Z1Z3 >. Equating the sum of 
the generators with a Hamiltonian, H = — J(ZiZ 2 + 
Z 2 Z 3 ), where J is the coupling strength, yields the Ising 
interaction between three spins in a line. If the full sta- 
bilizer is used, then H = -J{Z 1 Z 2 + Z 2 Z 3 + ZiZ 3 + I), 
which corresponds to Ising interactions between spins on 
a triangle with an energy offset due to the identity. In 
all cases, the codespace is the degenerate ground state of 
the Hamiltonian. 

We are interested in the thermodynamic limit and 
whether the energy gap between the ground state and the 
transition state will preserve the information. Our goal 
is to study how the choice of generators affects the crit- 
ical temperature of the ferromagnetic phase transition. 
At higher levels the choice of generators is non-trivial. 
As an example, consider k=2. The standard choice of 
generators is S k = 2 =< Z X Z 2 , Z 2 Z 3 , Z A Z b , Z b Z G , Z 7 Z 8 , 
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FIG. 1: Memory stabilizer structures generated by two body 
interactions. Circles are spin sites (qubits) and the lines 
show pairs of interacting spins (generators). The interaction 
strength J is constant (see Eq. [TJ. The total number of bits 
increases with concatenation level, k, as 3 fe . Only the first 
three levels are shown. 



ZgZg, Z\Z2Z 3 Z±Z^Z§, Z^Z^Z§ZiZ%Z§ >. Notice the six 
bit operators can be written as products of Z k=1 opera- 



k = l ryk = \ 



tors, Zf =1 Z. 



Z1Z2Z3Z4Z5ZQ. If the generators are 



used to define a Hamiltonian, then increasing k leads to 
many-body operators that act on 2 x 3 fe ~ 1 spins at once. 
This exponential increase in the many-body nature of 
the Hamiltonian makes the physical construction of such 
a system unlikely. 

In contrast, we can choose a set of pairwise Ising inter- 
actions that generate the same group. A natural choice 
would be the line S k = 2 =< Z X Z 2 , Z 2 Z 3 , Z 3 Z 4) Z4Z5, 
Z^Zq, Z5Z7, ZfZg, Z$Zg > but it is well-known that the 
ID Ising model does not have a phase transition at finite 
temperature. We instead consider a model where the en- 
coded bits are Ising coupled middle to middle not end 
to end: S k = 2 =< Z^ 2 , Z 2 Z 3 , Z 2 Z 5 , Z A Z 5 , Z 5 Z 6 , Z b Z 8 , 
Z 7 Z 8 , Z S Z 9 >. This describes the tree labeled Structure 
1 in Fig. [1] Structure 2 and 3 are modifications that 
include loops in the structure. The loops are equivalent 
to choosing a non-minimal set of generators. 

For each structure, the total number of bits, AT, in- 
creases with concatenation level, k, as N = 3 . The set 
of generators that corresponds to the structure defines 
the Hamiltonian, 



H — — J ^ ZiZj, 



(1) 



(id) 



where the sum is over nearest-neighbors. J sets the 
energy scale for the problem and temperature is mea- 
sured in units of J/ks- The stability of the structure 
as measured by the phase transition temperature, T c , 
depends on the energy-barrier that separates the two 
ground states and the number of pathways that traverse 
the barrier. For these complicated structures T c must be 
calculated numerically. 
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TABLE I: Average coordination and number of generator el- 
ements per spin site of the three memory stabilizer graphs 
described in Fig. [T]at the k — * oo limit. 



Structure 


Coordination number 


Generators per spin site 


1 


2 


1 


2 


3 


II 


3 



B. Calculating thermodynamic properties 

In the canonical ensemble, the thermody- 
namic average of an operator A, {A)t, is given 
by Tr[Aexp{-H/k B T)]/Z(T), where Z(T) 
Tr[exp(—H/ksT)] is the partition function. In 
practice this average cannot be exactly calculated in 
the limit of large N and numerical approximations are 
required. Our focus is on calculating the following 
thermodynamic properties: the average energy, (E)t, 
and the average absolute magnetization, (M)t- At T c 
the fluctuations of these quantities diverge for an infinite 
system, and T c can be determined by examining the 
specific heat capacity c(N, T) and the zero field magnetic 
susceptibility x(iV, T) as given by 



c(N,T) 



X(N,T) = 



1 <T 2 ) T - (E)* 



N 



k B T 2 



1 (M 2 ) T - (M) 2 T 
N k B T 



(2) 



(3) 



A direct measurement of the degree of preservation of the 
information can be read from the average magnetization 
per spin site, defined as m(N,T) = (l/N)(M)x- Below 
T c the system develops spontaneous magnetization and 
the single order parameter m approaches the value of 1. 

Without analytical expressions for (E)t and (M)t 
, there are three computational challenges: estimating 
thermodynamic averages for specific values of T, deter- 
mining T c from the evaluation of thermodynamic aver- 
ages at a finite set of T, and extending our results to the 
limit of large N. We solve each problem using well es- 
tablished numerical techniques for statistical mechanics. 

Reliable studies of thermodynamic averages near crit- 
ical temperatures require simulations of very large sys- 
tems. Two possible simulation methods are Metropo- 
lis Monte Carlo and the Wolff algorithm pH |2j|. A 
Metropolis Monte Carlo step updates the configuration 
of spins by flipping randomly (one at a time) N cho- 
sen spins. Groups of adjacent spins tend to point in the 
same direction near the critical region, giving rise to cor- 
relations in the system. The linear size of these clusters 
(correlation length, £) diverges at the critical tempera- 
ture and successive configurations of spins are generally 
strongly correlated. The efficiency of the Metropolis algo- 
rithm is hindered by the increasing of the number of steps 



needed to obtain uncorrelated spin configurations [22j]. 
One way to overcome this critical slowing down is by 
choosing the Wolff algorithm, in which a cluster is iden- 
tified and flipped at every Monte Carlo step. The size of 
the cluster is chosen to preserve detailed balance. The 
Wolff algorithm generates a Boltzman weighted set of 
spin configurations from where it is possible to calculate 
canonical thermodynamic averages. This cluster algo- 
rithm has previously been used to study systems with 
inhomogeneous local couplings such as the dilute Ising 
magnet [l8[ and Sierpinski carpets [2(| HH [H, HH and 
the efficiency of the Wolff algorithm seems to increase 
as the dimension is lowered [23j. As a result the Wolff 
algorithm is preferred over Metropolis for the structures 
in Fig.[TJ This is quantified in Section UlI CI 

In the limit of large N, T c will correspond to the tem- 
perature where the magnetic susceptibility is maximized 
T* ax . To calculate T* ax (N) for finite N, the histogram 
method is used. For a specific temperature, T>, the 
states randomly generated by the Wolff algorithm follow 
the Boltzmann distribution and can be used to calcu- 
late very good estimates of the thermal averages. The 
histogram method approximates the thermal averages at 
nearby temperatures by re-weighting the probability of 
choosing a spin configuration with exponential factors 
that account for the difference between the temperature 
of interest and Tq. The distance AT which can reliably 
be extrapolated away from Tq is given by [25| : 



AT 



T, 



1 



Ac(T )- 



(4) 



We find that it is safe to extrapolate ±2AT from the 
calculated central temperature, Tq. Going two standard 
deviations away from the mean sample energy still leaves 
5% of the samples in the region around the peak of the 
reweighted histogram. For a collection of a million inde- 
pendent spin configurations, 5% is 50000 samples which 
yield a reasonable estimate of the internal energy. 

Finally, the standard finite size scaling analysis devel- 
oped by Fisher [l3|, [25| is used to determine the critical 
exponents from the behavior of thermodynamic averages 
as a function of the system size measured in linear dimen- 
sion, L. According to the standard scaling hypothesis, 
and provided that the size of the system is large enough, 

the following scaling properties are expected at the Grit- 
cs — /3 7 

ical point: c oc L~ , m oc L~ and x L~ , where v is 
the correlation length exponent. The correlation length 
scales as £(T) oc where t = \T — T c \ /T c is a reduced 

temperature. 

For structures with well-defined dimension the linear 
size follows: 



L = N 1 ^, 



(5) 



The dimensions of the three structures of Fig. [T] are un- 
known. We assume that Eq. [5] holds and define v 1 = v ■ d 
as a correlation length exponent scaled to the system size. 
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Finite size effects replace the divergences at the critical 
point by finite peaks shifted away from T c . Effective crit- 
ical temperatures can thus be defined for each size and 
each physical quantity concerned (magnetic susceptibil- 
ity for example) as the positions of these maxima. The 
shift away from T c , to first order approximation can be 
written as 

T^^Tc + xo-N- 1 /"' (6) 

for the case of susceptibility. A fit of T£ ax against the 
system size N using Eqn. [6] gives an estimate of T c , xa 
and v' . 

Provided that T c and i/are known with a sufficient 
accuracy, the following power laws are observed at the 
critical point: c(N,T c ) oc , m(N,T c ) cx N~£ , and 
x(N,T c ) oc . The computation of the critical expo- 
nents a, /3, 7 can be deduced from the dependence on 
size of c, m and %, respectively. 

III. RESULTS AND DISCUSSION 

We have calculated T c for the three structures of Fig.[T] 
from finite-size scaling analysis of the magnetic suscepti- 
bility using concatenation levels k = 4 to k = 7. The 
thermodynamic averages were calculated from sets of 
millions of spin configurations generated by the Wolff 
algorithm. For each case, the magnetization autocor- 
relation function was calculated to find the number of 
successive cluster flips that separate independent spin- 
configurations. A fit of the auto-correlation function to 
an exponential reveals that for the structures studied here 
the autocorrelation time T ste p S is less than a single step 
ranging from 0.3 to 0.8 cluster flips. This is in contrast to 
the Metropolis method where initial attempts returned 
autocorrelation times between 300A and 5000iV possi- 
ble single-spin flips. Our thermal averages include every- 
other Wolff Monte-Carlo step after an initial thermaliza- 
tion period of 5 x 10 3 steps. Once T c is predicted, the 
size effects on c, m, x, and T s t ep s are studied. We find 
that the thermodynamics of finite structures can be de- 
scribed by critical scaling exponents and that the Wolff 
algorithm is efficient on these structures at the critical 
region. 

A. Finite size effects 

The positions of the effective temperatures, T^ ax (N) 
are first estimated by processing the data from short runs 
of 3 x 10 4 Wolff Monte Carlo steps from 0.05 to 2 ev- 
ery 0.05 T. A second set of short runs are performed 
over a region of 0.5 with 0.01 T resolution centered at 
guessed T* ax (N) values. Magnetic susceptibility maxi- 
mums and the corresponding effective temperatures are 
computed more precisely using the histogram method. 
Table [TT] shows the temperatures at which Wolff cluster 
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T [J/k B ] 

FIG. 2: Normalized magnetic susceptibilities as a function of 
T for different concatenation levels predicted using the his- 
togram method (Structure 3). The solid points are temper- 
atures at which Wolff cluster simulations are performed to 
obtain a set of 1 x 10 6 uncorrelated samples. Temperature 
ranges used in predictions are estimated from Eqn. [4] and 
shown in Table ITU 

simulations are performed to obtain a set of 1 x 10 6 un- 
correlated samples. The magnetic susceptibility is calcu- 
lated from Eqn. [3] and is re- weighted using the histogram 
method over the reliable temperature range estimated 
from Eqn. |U Figure [5] shows as an example the results of 
one of these experimental runs on Structure 3. We repeat 
this procedure five times for each level of concatenation 
and structure to check the reliability of the histogram 
method and to give error estimates on effective temper- 
atures (T* ax columns of Table [TT| . The T* ax monotoni- 
cally decreases with the system size for all structures. 

The values of T^ ax as a function of N are plotted 
in Fig. [3] for each structure. The solid lines are fits to 
the points using Eqn. [5] from which T c , xo an d v' are 
obtained (see Tabic IIIII) . Better estimates of these pa- 
rameters would require additional data points. Unfor- 
tunately, for k < 4, Eqn. [6] is no longer valid due to 
higher-order scaling corrections in the small N limit. Ad- 
ditional data points would require calculating the ther- 
modynamic properties at higher levels of concatenation. 
It may be possible to study bigger systems by using the 
Wang-Landau algorithm [26[ with a two dimensional en- 
ergy and magnetization joint density of states [27| . 

The results of the simulations show that one way to 
increase the critical temperature is by adding generators 
to each spin site. However, there is not a clear connec- 
tion between coordination number and T c (Tables UJ and 
as the critical temperature does not follow a 1, l|, 
lg progression when going from Structures 1 to 3. When 
adding generators to go from Structure 1 to Structure 
2 (looping only at the k = 1 level), the increase in T c 
is modest. Structure 3 has loops at all concatenation 
levels, and T c is almost doubled in comparison to the 
one of Structure 1. Another distinction of Structure 3 
is its higher symmetry. The extra bond changes the en- 
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TABLE II: Simulated temperature , T S i m , (confidence region 2 AT as per Eqn. [4| and the related T£ ax obtained from the 
Histogram method. The reported temperatures are in units of J jks and uncertainties quoted are 2a errors. 







T s im 








T x (N) 

' max v- 1 * / 




Structure 


k = 4 


k = 5 


k = 6 


k = 7 


k = 4 


k = 5 k = 6 


k = 7 


1 


0.735(251) 


0.685(137) 


0.645(76) 


0.620(43) 


0.736(1) 


0.687(1) 0.650(5) 


0.611(6) 


2 


0.785(322) 


0.725(182) 


0.675(105) 


0.640(60) 


0.783(1) 


0.723(2) 0.678(3) 


0.646(3) 


3 


1.240(394) 


1.155(216) 


1.105(122) 


1.075(70) 


1.239(1) 


1.162(4) 1.099(6) 


1.056(8) 



-0.02 




FIG. 3: T^ a;c as a function of iV for each of the memory 
stabilizers of Fig [T] The solid lines are fits to the points 
using Eqn. [6] (fitting parameters T c , xo an d v are reported in 
Table Uri}. 



FIG. 4: Magnetizations at critical temperature against sys- 
tem size N for each memory stabilizer on log- log scale. Solid 
lines show least square fits to power laws from where the /3/v ' 
exponents are calculated. 



TABLE III: Finite size scaling law results (T c , xo, v') using 
data from Table [TT] in Eqn. [6] for the three structures shown 
in Fig. [T] Uncertainties quoted are 2a errors. 

Structure T c xo 1/1 

T 0.455(111) 0.603(35) 5.747(2.610) 

2 0.552(16) 0.769(41) 3.648(367) 

3 0.890(73) 0.953(69) 4.374(1.199) 



TABLE IV: Exponents obtained from power law behavior of 
magnetization and susceptibility at critical temperature for 
each memory stabilizer structure. Uncertainties quoted are 
2a errors. 



Structure 


0/u' 


7/1/ 


1 


0.015 


0.967(7) 


2 


0.032 


0.944(2) 


3 


0.027 


0.952(2) 



ergy gap between differing spin configurations and, due 
to symmetry, changes the underlying density of states at 
a given energy. 

We attempted to study the finite size effects on the 
heat capacity to further validate T c predictions. A wide 
fluctuation of the position of T^ ax from experiment to 
experiment was observed. It is hard to follow trends from 
the heights of specific heat peaks for different system sizes 
and structures. Bhanot et al. (28j pointed out that when 
space dimensionality is lower than 2, a is expected to be 
negative, and the specific heat versus temperature peak 
broadens as the system size increases. We were unable 
to extract ajv' in a reliable way from fits using T^ lax 
and c. Monceau and Perreau encounter similar problems 
on fractal structures of dimensionality between one and 
two [2l|. 



B. Magnetization and magnetic susceptibility at T c 

The mean values of the magnetization and the zero- 
field susceptibility are obtained from simulations at the 
previously calculated critical temperatures T c (Tablo lDl)) . 
We use a single set of 1 x 10 6 uncorrelated samples for 
the analysis in this section. The power laws m(N, T c ) oc 

and x( N , T c) oc are satisfied. Figures H and M 
show plots of the average absolute magnetization per spin 
and magnetic susceptibility versus system size TV on a 
log-log scale respectively. The results of the least square 
analysis are displayed in Table IIVI 

The exponent a is deduced from the Rushbrooke scal- 
ing law a — 2 — 2/3 — 7 an d it takes, as expected, negative 
values for the three memory stabilizers. Table [V] shows 
the set of critical exponents (cv, (3, 7) for each struc- 
ture. We write the Rushbrooke and Josephson scaling 
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FIG. 5: Magnetic susceptibilities at critical temperature 
against N for each stabilizer structure on log-log scale. Solid 
lines show least square fits from where the 7/z/ exponents are 
calculated. 



TABLE V: Set of critical exponents (a, j3, 7) of three struc- 
tures shown in Fig. [T] Last column checks consistency of the 
results by using the Rushbrooke and Josephson's scaling laws 
as discussed in the text. Uncertainties quoted are 2<r errors. 



Structure 



n 



7 



I = —— 



2l 



1 -3.730(2.526) 0.086(39) 5.558(2.525) 0.997(3) 

2 -1.679(347) 0.117(12) 3.444(347) 1.009(1) 

3 -2.395(1.143) 0.116(32) 4.163(1.141) 1.005(1) 



law d = j/v + 2f3/v as a function of v' to get rid of 
the unknown dimension d. The last column of Table IVl 
shows that the scaling law is satisfied within an error of 
less than 1%. This means that the magnetization (the 
order parameter of interest for memory preservation) is 
continuous at T c . The uncertainty in the critical expo- 
nents is quite big for the three graphs. It is not possible 
to conlcude whether or not the memory structures share 
the same set of critical exponents. All three structures 
may be in some weak universality class in which critical 
exponents may not only depend upon the symmetry of 
order parameters and fractal dimensions, but also upon 
their geometry. This seems to be the case for Ising mag- 
nets from Sierpinski fractals of non-integer dimensions 
between one and three Sol, EL HI fH . 



C. Wolff algorithm efficiency at T c 

The dynamical aspects of the Wolff algorithm when 
applied to memory stabilizers are analyzed. We take five 
runs of 5 x 10 5 cluster flips at each concatenation level 
k = 4 — 7 and for each structure to calculate the magneti- 
zation autocorrelation function. The autocorrelations are 
fit to a single exponential decay to obtain Wolff autocor- 
relation times T s teps (see Table I VI I) . As shown in Fig. [51 
magnetization autocorrelation times follow the power law 
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FIG. 6: Magnetization Wolff autocorrelation times (in num- 
ber of cluster flips) against system size N for each memory 
stabilizer structure on log-log scale. The solid lines are fits to 
T s t e ps oc N z °/ d from where measured critical exponents zo/d 
are calculated. 



TABLE VI: Wolff autocorrelation times, T s t eps , measured 
from the decay of the magnetization autocorrelation function 
for each memory stabilizer structure and system size. Last 
two rows show the measured critical exponent zq /d, obtained 
from power law fits, and the Wolff dynamical critical exponent 
z/d (from Eqn. [7J|- Uncertainties quoted are 2<r errors. 



k 


1 


Structure 
2 


3 


4 


0.284(16) 


0.296(26) 


0.413(18) 


5 


0.326(12) 


0.366(20) 


0.465(19) 


6 


0.359(11) 


0.434(18) 


0.508(16) 


7 


0.390(10) 


0.492(18) 


0.557(18) 


zo/d 


0.090(15) 


0.145(26) 


0.088(8) 


z/d 


0.057(16) 


0.089(26) 


0.040(8) 



T steps oc N z °/ d at the critical temperature. The Wolff dy- 
namical critical exponent, z/d, associated with memory 
stabilizer structures is defined as 



4 = ^ + ^-1. 



(7) 



The Wolff algorithm is very efficient in reducing the 
critical slowing down (increase in correlation time as T c 
is approached) for the stabilizer structures of Fig. [T] The 
dynamical critical exponents are very low compared to 
the Metropolis or Wolff algorithm on the 2D Ising Model, 
where z/d= 1.0835 and z/d — 0.125 respectively [25| . 



IV. CONCLUSIONS 

We have used the Wolff algorithm together with the 
histogram method and finite-size scaling analysis to cal- 
culate critical temperatures of Hamiltonians based on 
concatenated error-correction codes. The three simple 
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two-body-interaction structures investigated have differ- 
ent leveis of connectivity. We find that the relationship 
between coordination number and critical temperature 
is not obvious. The intriguing result is that the number 
of generators is less important than the structure. For 
a minimum number of generators, one can have either a 
linear Ising model with no phase transition or Structure 
1 with a finite phase transition. If one adds additional 
connections or generators, the results can range from a 
modest increase in T c (Structure 1 and Structure 2) to a 
doubling of T c (Structure 1 and Structure 3). Whether 
these insights can be applied to self-correcting quantum 
systems is an open question. 

Scaling properties of the magnetization and magnetic 
susceptibility satisfy power-law fits as function of total 
number of spins N. Each structure exhibits second or- 
der or continuous phase transition. We report the set of 
critical exponents (a, /3, 7), and by fitting the decay of 
the magnetization autocorrelation functions at the crit- 
ical points we calculate Wolff dynamical critical expo- 
nents. It is possible that all three structures are in some 
weak universality class but the current study does not 



show this. 

For quantum information, a thermodynamically unsta- 
ble memory that has a large kinetic barrier could be use- 
ful for preserving information. It is possible that many 
of the stabilizer codes that are not self-correcting mem- 
ories could satisfy this relaxed condition. Although the 
kinetics depends strongly on the details of the specific 
system-bath coupling [29j |. the work here suggests that 
the choice of geometry and generators could lead to large 
differences in the effective information preservation. 
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